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Abstract 

Using data collected near the D* + D~ peak production energy E cm = 4170 MeV by the CLEO-c 
detector, we study the decays of Df mesons to two pseudoscalar mesons. We report on searches 
for the singly-Cabibbo-suppressed D + decay modes K + rj, K + r]', ir + Kg, K + ir°, and the isospin- 
forbidden decay mode L> + — ► 7r + 7r . We normalize with respect to the Cabibbo-favored modes 
7r + 7/, 7T + rj', and K + Kg, and obtain ratios of branching fractions: B(Df — ► K + 7]) / B(Df — > 7r + r/) 
= (8.9 ± 1.5 ± 0.4)%, B(D+ -► K+r?') / -► vr+r/) = (4.2 ± 1.3 ± 0.3)%, B(D+ -► vr + K°) 

/ -»• K+Kg) = (8.2 ± 0.9 ± 0.2)%, £(£>+ -► K+ir°) / B{Df -► K+Kg) = (5.0 ± 1.2 ± 

0.6)%, and fi(L>+ -► vr+vr ) / B(D+ -► K + K^) < 4.1% at 90% CL, where the uncertainties are 
statistical and systematic, respectively. 
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There are ten possible decays of Df mesons to a pair of mesons from the lowest-lying 
pseudoscalar meson nonet. The decay can be to either K + or 7r + , combined with any of 77, 
7/, 7T°, K°, or K° (K s or K\ for the final state). Measurements of the branching fractions 
of the complete set of decays test flavor topology and SU(3) predictions [lj]. The Cabibbo- 
favored, color-favored (external spectator) decays — > 7i + r] and — > 7t + t/ have been 
previously measured |2j, as has the Cabibbo-favored, color- mixed (internal spectator) decay 
Df — > K + K s [2]. Here we present first observations of the singly-Cabibbo-suppressed, 
color-favored decays Df —>■ K + rj, — > K + rj', and Df — > n + K s , and strong evidence (4.7 
standard deviations (a)) for the singly-Cabibbo-suppressed, color-mixed decay K + ir°. 
(In this analysis, we have detected K s , but made no attempt to detect K®, nor have previous 
Df measurements.) We measure the ratio of the branching fraction of each singly-Cabibbo- 
suppressed decay to that of the corresponding favored decay, expected to be, and found to 
be, of order |V^/V^ S | 2 ~ 1/20. The decay — ► 7r + 7r° requires a change in isospin of 2 
units, and is thus "isospin- forbidden" , and expected to be substantially suppressed. Our 
search for this decay reveals no firm evidence for it, and we present an upper limit. 

CLEO-c is a general-purpose solenoidal detector. The charged particle tracking system 
covers a solid angle of 93% of 4tc and consists of a small-radius, six-layer, low-mass, stereo 
wire drift chamber, concentric with, and surrounded by, a 47-layer cylindrical central drift 
chamber. The chambers operate in a 1.0 T magnetic field and achieve a momentum res- 
olution of ~0.6% at p =1 GeV/c. We utilize two particle identification (PID) devices to 
separate charged kaons from pions: the central drift chamber, which provides measurements 
of ionization energy loss (dE/dx), and, surrounding this drift chamber, a cylindrical ring- 
imaging Cherenkov (RICH) detector, whose active solid angle is 80% of 47r. Detection of 
neutral pions and eta mesons relies on an electromagnetic calorimeter consisting of 7784 
cesium iodide crystals and covering 95% of 4ir. The calorimeter achieves a photon energy 
resolution of 2.2% at £ 7 =1 GeV and 6% at 100 MeV. The CLEO-c detector is described 
in detail elsewhere jjjj]. 

We use 298 pb _1 of data produced in e + e~ collisions using the Cornell Electron Storage 
Ring (CESR) near the center-of-mass energy a/s = 4170 MeV. Here the cross-section for 
the channel of interest, D* + D~ or D~£~D*~, is ~1 nb [4j]. We select events in which the D* 
decays to D s + 7 (94% branching fraction [2J). Other charm production totals ~7 nb (H, 
and the underlying light-quark "continuum" is about 12 nb. We reconstruct mesons in 
all two-body pseudoscalar decay channels. Throughout this Letter, charge conjugate modes 
are implicitly assumed, unless otherwise noted. 

We use the reconstructed invariant mass of the D s candidate, M(D S ), and the mass 
recoiling against the D s candidate, M Tecoi \(D s ) = J (y/s — E Ds ) 2 — pf^, as our primary kine- 
matic variables to select a D s candidate. Here f>o a is the momentum of the D s candidate, 
Ed s = J m 2 Ds + p 2 Da , and mn s is the known D s mass [2]. We make no requirements on the 
decay of the other D s in the event. 

There are two components in the recoil mass distribution, a peak around the D* mass 
if the candidate is due to the primary D s and a rectangular shaped distribution if the 
candidate is due to the secondary D s from D* decays. The edges of M Teco n(D s ) from the 
secondary D s are kinematically determined (as a function of y/s and known masses), and 
at v/i = 4170 MeV, AM recoil (D s ) = M recoil (D s ) - m D * is in the range [-54, 57] MeV. Initial 
state radiation causes a tail on the high side, above 57 MeV. We select D s candidates within 
the -55 MeV < AM recoil (Z) s ) < +55 MeV range. 

We also require a photon consistent with coming from D* + — > Dfj decay, by look- 
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ing at the mass recoiling against the D s candidate plus 7 system, M iecoi i(D s + 7) = 

\J (a/s — E Ds — E^) 2 — (pd s + P-y) 2 - For correct combinations, this recoil mass peaks at mc,, 
regardless of whether the candidate is due to a primary or a secondary D s . We require 
\M reco \\(D s + 7) — moj < 20 MeV. Though there is a 25% efficiency loss from this require- 
ment, it improves the signal to noise ratio, important for the suppressed modes. 

Our standard final-state particle selection requirements are described in detail else- 
where jEf. Charged tracks produced in the decay are required to satisfy criteria based 
on the track fit quality, have momenta above 50 MeV/c, and angles with respect to the 
beam line, 9, satisfying | cos^l < 0.93. They must also be consistent with coming from the 
interaction point in three dimensions. Pion and kaon candidates are required to have dE/dx 
measurements within three standard deviations (3<r) of the expected value. For tracks with 
momenta greater than 700 MeV/c, RICH information, if available, is combined with dE/dx. 
The efficiencies (95% or higher) and misidentification rates (a few per cent) are determined 
with charged pions and kaons from hadronic D decays. 

The Kg candidates are selected from pairs of oppositely-charged and vertex-constrained 
tracks having invariant mass within 12 MeV, or roughly 4.5a, of the known Kg mass. We 
identify tt° candidates via tt° — > 77, detecting the photons in the Csl calorimeter. To avoid 
having both photons in a region of poorer energy resolution, we require that at least one of 
the photons be in the "good barrel" region, | cos# 7 | < 0.8. We require that the calorimeter 
clusters have a measured energy above 30 MeV, have a lateral distribution consistent with 
that from photons, and not be matched to any charged track. The invariant mass of the 
photon pair is required to be within 3a (a ~ 6 MeV) of the known 7r° mass. A 7r° mass 
constraint is imposed when 7r° candidates are used in further reconstruction. We reconstruct 
77 candidates in two decay modes. For the decay rj — > 77, candidates are formed using a 
similar procedure as for 7r° except that a ~ 12 MeV. For 77 — > 7r + 7r~7r°, we require that the 
invariant mass of the three pions be within 10 MeV of the known 77 mass. We reconstruct 
rj candidates in the decay mode rj — > 7r + 7r~77. We require \m n + n - v — m v '\ < 10 MeV. 

The D s invariant mass distributions of the backgrounds to £)+ — > K + K s and — > 
n + K s are not smooth, but have bumps, caused by D* + D*~ events followed by D* ± — > n^D 
decays. The low-momentum tc^ from D*^ decay, in combination with a particle from D° 
decay, can create a fake Kg. To reduce the bump structure, which complicates fitting the 
background, we reject those — > K + K s and — > tt + K s candidates that contain a n + or 
7r~ with momentum below 100 MeV/c. For symmetry, we also reject events with a K ± with 
momentum below 100 MeV/c. Further, we require that the Kg has traveled a measurable 
distance from the interaction point before decaying, i.e., that the distance along the flight 
path, from interaction point to K s decay vertex, be greater than zero with a 3a significance. 
After the low-momentum track veto and K s flight significance requirement are applied, no 
bump structures remain. 

For the modes with rj or rj', 77 — > 77, we reject the r] candidate if either of the daughter 
photons is consistent with coming from 7r° — > 77 when paired with any other 7 in the event. 
This veto reduces the background from fake 77 combination for rj candidates. 

The resulting M(D S ) distributions for the Cabibbo-favored and Cabibbo-suppressed D s 
modes are shown in Figs. [1] and [21 respectively. The points show the data and the lines 
are fits. We perform a binned maximum likelihood fit (2 MeV bins) to extract signal yields 
from the M(D S ) distributions. For the signal, we use the sum of two Gaussians for the 
line shape. The signal shape parameters are determined by fits to M(D S ) distributions 
obtained from a GEANT-based Monte Carlo (MC) simulation 0. For the background, we 
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FIG. 1: M(D S ) distributions for Cabibbo-favored D s modes from data. The points are the data 
and the superimposed line is the fit (the dotted line is the fitted background) as described in the 
text. 

use a second-degree polynomial function, allowing the overall scale, and the coefficient of 
the linear term relative to the constant term, to float in the fits to the data. We constrain 
the (very small) coefficient of the quadratic term relative to the constant term to the value 
given by MC simulation. We include as a systematic error the change in yield caused by 
varying the quadratic coefficient over a reasonable range, typically doubling the quadratic 
term coefficient, or setting it to zero. (For the favored modes, where the background is 
relatively smaller, we allow the coefficient of the quadratic term to float.) 

Results of the fits are shown in Table [B Also given in Table [I] is the detection efficiency 
for each mode, and, for the Cabibbo-suppressed modes, the statistical significance of the 
signal. For — > K + tt°, the statistical significance is 4.7 standard deviations (a), while for 
the other modes using rj — > 77, the statistical significance exceeds 5a. The 77 — > 7r + 7r~7r° 
mode for — > K + t] confirms the signal, at 3.7a, while for Df — > K + i]', due to the very 
large background of this mode, it gives no supporting evidence. For the D+ — > ir + Kg mode, 
the statistical significance exceeds 10a. For all Cabibbo-favored modes, very clear signals 
are found in the data. 

We find no significant evidence for the isospin-forbidden decay — > 7r + 7r°, and therefore 
set an upper limit on its rate. There is a large background from continuum events, and Monte 
Carlo studies indicate that tightening the requirement on recoil mass to ±10 MeV should 
improve the upper limit. The invariant mass distribution with this requirement applied 
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FIG. 2: M(D S ) distributions for Cabibbo-suppressed D s modes from data. The points are the 
data and the superimposed line is the fit (the dotted line is the fitted background) as described in 
the text. Also shown is the distribution for the isospin-forbidden decay Df — > 7r + 7r°. 

is shown in Fig. [2j We apply a sideband subtraction to the invariant mass distribution 
and obtain a yield of 17 ± 25 events. We interpret this result as implying a probability 
distribution for the true number of events N as a Gaussian, centered on 17, with width 
a = 25, but truncated at zero, so the probability distribution vanishes for a negative true 
number of events. Ninety percent of the area of this distribution lies below 52 events, 
which we take as the 90% confidence level upper limit on the true number of events. We 
normalize this upper limit on yield to that for — > K + K s , obtaining B(D+ — > 7r + 7r°) 
/ B(Df -> K+K° s ) < 3.80 x 10~ 2 (statistical only). Systematic errors, from the ratio of 
detection efficiencies, are ±1.8% for the Kg, ±4.2% for the tt , and other smaller errors, 
leading to a combined relative systematic error of ±5.1%. We conservatively increase the 
upper limit by 1.28 times the combined systematic errors, giving an upper limit, including 
systematic errors, of 4.1 x 10~ 2 . 
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TABLE I: Observed yields from data and reconstruction efficiencies and their statistical uncer- 
tainties. For the Cabibbo-suppressed modes, the statistical significance of the signal is also given 
(see the text for details). The efficiencies include sub-mode branching fractions [2j, and have been 
corrected to include several known small differences between data and Monte Carlo simulation. 
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4.39 ± 0.03 


Df 


K + rf 


rf- 
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77 


25 ± 7 


5.7 


2.10 ± 0.02 


Df 




rf- 


* ir + ir~r], rj — ► 


7T + 7T _ 7r 


3 ± 6 


0.0 


1.53 ± 0.01 



TABLE II: Ratios of branching fractions of Cabibbo-suppressed modes to corresponding Cabibbo- 
favored modes. Uncertainties are statistical and systematic, respectively. 
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Mode B S /B F (10 

B{Df -»■ K+r]) I B{Df -»■ vr+r/) 8.9 ± 1.5 ± 0.4 

B{Df -» K+rf) I B(Df -» TT+r?') 4.2 ± 1.3 ± 0.3 
0(D+ -► 7r+K§) / S(I>+ -> A+A£) 8.2 ± 0.9 ± 0.2 
£(£>+ -» A+7T ) / B(Df -> A+A°) 5.0 ± 1.2 ± 0.6 
£(£>+ -► tt+tt ) / B(£>+ -> A+A°) < 4.1 (90% CL) 



In principle, non-resonant D s decay could enter into our signal modes with the same 
final particles. For example, non-resonant Df — ► 7r + (7r + 7r~7r°) could appear in the Df — > 
7i + i], 7] — > 7T + 7r _ 7r° mode. To understand the background from non-resonant D s decay, we 
look at M(D S ) distributions in the sideband region of the intermediate resonance (77, rj', or 
Kg) invariant mass. These studies show that the non-resonant modes produce negligible 
contributions to our signal modes. 

For the modes with 77 or rj' (77' — > 7T + 7T~r]) in the final state, we reconstruct these modes 
with 77 decaying to 77 and to 7r + 7r~7r°. For Cabibbo-favored modes, we combine the two 
fit yields from the different T] decay modes according to the fit yield fractional error. The 
weighting factors for both Df — > tt + i] and Df — > 7r + r/' are 0.65 for 77 — > 77 and 0.35 
for 77 — > 7T + 7r _ 7r . We apply the same weighting factors to the corresponding Cabibbo- 
suppressed modes (Df — > K + r] and Df — > K + r)'). Doing so guarantees cancellation of 
systematic errors between Cabibbo-favored and Cabibbo-suppressed modes. It also avoids 
a possible bias that could come from using the errors on the Cabibbo-suppressed modes to 
determine the weighting factors for them. 

Ratios of branching fractions are computed for each of the Cabibbo-suppressed modes 
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TABLE III: Measured CP asymmetries in Cabibbo-suppressed decay modes. Only statistical 
uncertainties are included. Systematic errors are negligible by comparison. 



Mode 




(B + -B_)/(B + +B_)(%) 


A{D+ - 


-» K+rf) 


-20 ± 18 


A(D+ - 


-> K+rf) 


-17 ± 37 


A(Df - 


-> vr+K°) 


27 ± 11 


A(D+ - 


-> K+7T ) 


2 ± 29 



and are presented in Table [TTJ They are normalized with respect to the corresponding 
Cabibbo-favored modes. We use the —>■ K + K s mode to normalize the — ► K + tt° 
mode. The upper limit for the unobserved mode — ► 7r + 7r , normalized with respect to 
— > is also shown in Table [TTl 

We have considered several sources of systematic uncertainty. Finite MC statistics in 
determining reconstruction efficiencies introduces uncertainties at the level of less than 1%. 
The uncertainty associated with the efficiency for finding a track is 0.3%; an additional 0.6% 
systematic uncertainty for each kaon track is added. The relative systematic uncertainties 
for 7T° and K s efficiencies are 4.2% and 1.8%, respectively. The systematic uncertainty for 
7] efficiencies cancels in all ratios. Uncertainties in the charged pion and kaon identification 
efficiencies are 0.3% per pion and 1.3% per kaon j^]. The systematic uncertainties from 
the Kg flight significance requirement and the low-momentum track veto are 0.5% and 
0.3%, respectively. The signal shape parameters are taken from MC simulation, and have 
uncertainties related to possible flaws in simulation. We estimate this systematic uncertainty 
by allowing the signal shape parameters of the favored modes to float, and having the 
parameters of the suppressed modes track those of the favored modes. We find uncertainties 
of 0.2% to 2.9%, depending on mode. For the suppressed modes, the background quadratic 
term is also taken from MC simulation. We vary that term over a reasonable range, finding 
a systematic error of 2.4% to 10.6%, depending on mode. 

In calculating the relative systematic uncertainties for the measured ratio of Cabibbo- 
suppressed mode branching fractions to Cabibbo-favored mode branching fractions 
( ^Suppressed /^Favored), cancellation of uncertainties has been taken into account. The sys- 
tematic uncertainties that do not cancel in the ratios are added in quadrature to obtain 
the total systematic uncertainties shown as the second error in Table [TTl For the upper 
limit in Table [TTJ the systematic uncertainties have been included as previously described. 
Systematic uncertainties for all measured ratios are at most half the statistical uncertainties. 

The Standard Model predicts that direct CP violation in D decays, e.g., a difference 
in the branching fractions for Z?+ — > K + r) and Dj — > K~t], will be vanishingly small. As 
a search for evidence of non-Standard-Model physics, we have therefore measured the CP 
asymmetries A = (£> + — £?_)/ (£> + + £L) for the four Cabibbo-suppressed D s decay modes we 
are studying. Results are given in Table IIHI Errors shown are statistical. The systematic 
errors, from the differences in efficiency for detecting K + vs. K~ and tt + vs. tt~ , are < 2.0%, 
negligible by comparison. All asymmetries are consistent with zero. 

In summary, we report first observations of four Cabibbo-suppressed decays of D s mesons, 
and measure the ratio of their branching fractions to the corresponding Cabibbo-favored 
modes. We find those ratios to be of order IV^/V^I 2 ~ 1/20 in agreement with naive 
expectations. We report a first upper limit on the isospin-forbidden decay — > 7r + 7r°. 
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The CP asymmetries for the four Cabibbo-suppressed decays are consistent with zero, as 
predicted by the Standard Model. 
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